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We numerically analyze spectral properties of the Fibonacci model which is a one-dimensional 
quasiperiodic system. We find that the energy levels of this model have the distribution of the 
band widths w obeys Pb{w) ~ w a (w — * 0) and Pb{w) ~ e~ l3w (w — > oo), the gap distribution 
Pa{s) ~ s~ s (s -> 0) (a, (3,5 > 0) . We also compare the results with those of multi-scale Cantor 
sets. We find qualitative differences between the spectra of the Fibonacci model and the multi-scale 
Cantor sets. 
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I. INTRODUCTION 

Bloch's theorem tells that the free electron eigen- 
states are always extended in periodic systems. It is also 
well known that the free electron eigenstates are always 
localized in one-dimensional random media 0. On the 
other hand, one-dimensional quasiperiodic systems show 
various interesting behavior 0] • A one-dimensional tight- 
binding model is 

ti+li>i+l + U-iipi-i + ZiVJi = Etpi, (1) 

where ipi denotes the value of the wave function at i- 
th site, ti and ei are the hopping matrix element and 
the site energy at i-th site respectively, both of which 
can be taken to be quasiperiodic with lattice spacing. 
The Harper model: ti — 1 and e< = Acos(27reri + 6) 
has been well studied. When a is an irrational number, 
it is quasiperiodic. The eigenstates of the quasiperiodic 
Harper model are extended for A < 2 and are localized 
for A > 2 with the metal-insulator transition point A = 2 
@ . Hofstadter found a surprisingly rich structure of the 
spectrum at A = 2 pj . There the spectrum as well as the 
eigenstates become multifractal |2j. The total measure 
of the bands at the critical point is zero with a fractal 
dimension less than one 5,6]. The scaling behavior of the 
spectrum has been extensively studied |2| by multifractal 
analysis 0. 

Recently a new statistical characterization of the spec- 
trum of the Harper model at the critical point was found 
0, || . The distribution of the normalized width of the 
bands was examined. The distribution of the band- 
widths Ps(w)dw counts the number of bands whose val- 
ues are between w and w + dw. In [8|, Pb(w) was nu- 
merically investigated at A = 2 and a — v ^~ 1 and it was 
claimed that Pb(w) follows the semi-Poissonian distribu- 
tion 4we~ 2w . However this claim has been excluded by 
further detailed numerical investigations Q. Still it has 
been confirmed that 

P B (w)~w a (W-.0), (2) 



and 

Pb{w) ~ e- pw (w -> oo), (3) 

where a ~ 2.5,0 ~ 1.4 |. Similar laws are also found 
for a variant of the Harper model at its criticality Q . 

For the Harper model, the distribution of energy gaps 
Pg(s) was also examined. The distribution div erge s near 
the origin and follows an inverse power law 0, Il0| 

P G (s)~s- s (s-0), (4) 

with 5 ~ 1.5. This law is also extended to a variant of 
the Harper model 

In order to investigate how these laws are univer- 
sal for critical quasiperiodic systems, we study the one- 
dimensional Fibonacci model [ll], [lj ■ It is a model ob- 
tained by setting ti or ei in (JIJ to the Fibonacci sequence. 
We shall study the case with = 0. The Fibonacci 
model remains critical when varying the parameters ti, 
and the spectra are multifractal similar to that at the 
critical point A = 2 of the Harper model Q . We will find 
qualitatively similar but quantitatively different behavior 
oiP B (w) and P G {s)- 

A remarkable feature of this model is that the con- 
struction of the Fibonacci sequence can be translated 
to a renormalization group transformation, and its ac- 
tion on the trace of the transfer matrix becomes a non- 
linear dynamical system called the Kohmoto-Kadanoff- 
Tang (KKT) map[n|. The KKT maphas been analyzed 
as a nonlinear dynamical system [ljl llJ, [lj| . Namely the 
periodic orbit corresponding to the bands' center was an- 
alyzed, and the scaling property of the band was shown 
to be determined by the eigenvalues of the linearized map 
at the hyperbolic points. They induce heteroclinic points 
of the KKT map. The structure of subsets of the spec- 
trum was explained by the Smale's horseshoe structure. 
This analysis can be extended to other bands and then 
the explanation is generalized to a subset of the whole 
spectrum [l5| . 

The Smale's horseshoe structure indicates that the 
spectrum of the Fibonacci model is a multi-scale Cantor 
set. Here, we restrict the meaning of the term "multi- 
scale Cantor set" as the set which is obtained by the 
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following procedure : we transform the interval I = [0, 1] 
into a finite number of subintervals with definite scales 
and apply the transformation to subintervals successively. 
In the literature, the term "Cantor set" is also used as 
the term for "a closed set with no isolated points and 
whose complement is dense". We call this "general Can- 
tor set" . We will compare our results with those of multi- 
scale Cantor sets. Also note that the spectrum of the 
Fibonacci model is proved to be singular continuous 0] . 
Nevertheless the comparison will turn out to be useful. 

In Sec. II, we recall the renormalization group trans- 
formation. In Sec. Ill, we report the results of our nu- 
merical studies. In Sec. IV, wc compare our results with 
those of multi-scale Cantor sets. Sec.V gives conclusion. 
In Appendix, we provide definition and distributions of 
multi-scale Cantor sets. 



II. FIBONACCI MODEL AND THE 
KOHMOTO-KADANOFF-TANG MAP 
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We recall the aspects of the Fibonacci model we use, 
especially the renormalization group transformation . 

The Fibonacci model describes a one-dimensional 
quasiperiodic system given by JQ) with e.; — 0, defined 
by 

ti+itpi+i + tiipi-i = Etpi, (5) 

where t^s take two values t a and % arranged in the 
Fibonacci sequence constructed recursively as S^+i = 
{Sg,Sg—±} (£ > 1) with an initial condition So = 
{b},Si = {a}. The number of elements in Sg is the Fi- 
bonacci number Fg defined by Fg+i = Fg + Fg_i (£ > 1) 
with an initial condition Fq = F± = 1. To investigate 
the model with infinite sites, one considers a finite size 
system of © with Fg sites and imposes a boundary con- 
dition on the states. Then the system is equivalent to 
the system with a period Fg. 
We write Eq. © as 

(t'Mf'r )(£)■ (6) 

The present 2x2 matrix is a transfer matrix. Values 
of the wave function at any lattice sites are related to 
those of the initial sites by applying the transfer matrices 
successively. There are three types of transfer matrices : 

M " = (f "o 1 )' 

/ E_ _tt_ \ / E _ta \ 

Mab = \ J- J , M ba = h h \ . (7) 

We introduce a matrix Mg which generates a wave func- 
tion at the Fg-th site 

(*&')-*(£)• < 8 » 



FIG. 1: Two-dimensional manifold x 2 +y 2 + z 2 — 2xyz — 1 = / 
for I — The manifold consists of a central body and four 
parts to infinity. The four parts are connected to the central 
body through four tubes. 



This matrix Mi is given by a successive product of the 
transfer matrices, and the number of element of the prod- 
uct is Fg which grows like a~ l . A special property of the 
transfer matrices of the Fibonacci model is that Mg is 
determined by a recursive map pd| 

M e+1 = Mg-tMg, (9) 

with an initial condition M\ — M aa , — M^M^a. This 
map simplifies a procedure to obtain the matrix Mg. 

As Mg's are unimodular matrices, using the trace of 
the matrix Mg, the equation © is casted to 

xg+i = ye, ye+i = zg, z i+1 = 2ygzg - xg, (10) 

where xg = ^Tr Mg and we introduce three-dimensional 
vectors fg = (xg,yg,zg) = (xg,xg-i,xg- 2 ) .11]. This map 
Fe+i = fi^t) is the KKT map. An initial condition is 

xi = ^f-,?/i = j^i z i = \ {^t~ + ff)- ^ i s remarkable 
that this map has a constant of motion |13| : 

I = xj+yj + zj - 2x e y e zg - 1. (11) 

The initial condition gives the value of this constant of 
motion as 

I= l( t A^ t A\ (12) 
4 \t a t b J K ' 

The constant of motion / determines a two-dimensional 
manifold on which an orbit of the KKT map remains. In 
Fig. 1, an example of the manifold is shown. It is al- 
ways non-compact as long as / > 0. When / approaches 
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FIG. 2: Energy bands at (t a ,t b ) = (1,2), I = ^ with Fi- 
bonacci numbers F3 = 3, F4 — 5, F5 = 8 and Fe = 13. The 
bands consist of Fi bands and F e — 1 gaps. 



FIG. 3: Band-width distributions at {t a ,h) = (1,2), I = ^. 
The distributions for lniu < 1.5 suggest Pb(w) ~ iu 1 ' 71 . 



zero, the tubes shrink to points. A case / < does not 
represent the electronic systems. 

In Fig. 2, we show the energy bands for small £'s. 
Since the eigenstates of the model must be normalizable, 
we impose the condition on the initial condition so that 
\ipFi\ does not grow exponentially. In other words, the 
energy levels of the equation JSJ with Fg is determined 
by the condition \xg\ < 1 since the energy enters in the 
initial condition of the KKT map. It can also be shown 
that a condition \xi\, > 1 for an i < £ is sufficient 

for the corresponding energy to be in an energy gap. For 
large Fe, the KKT map provides escaping orbits on the 
manifold for almost all energies and the remaining set is 
shrinking to zero measure as Fe — * 00. 



III. LEVEL STATISTICS 



„ -1.5 
5 



I 


i i i 

17711 


1 

□ 




28657 






46368 







i i i 


1 



0.5 1 1.5 2 2.5 3 3.5 
w 

FIG. 4: Band- width distributions with w-1iiPb(w) plots at 
(t a ,tb) = (1, 2), I = jjj. The distributions for w > 0.5 suggest 
P B (w) ~ e - 1M8w . 



We use the KKT map described in the previous sec- 
tion to get the spectrum of the Fibonacci model. We first 
consider periodic systems whose unit cell is a Fibonacci 
number Fe, then extrapolate results to understand the 
Fibonacci model £ = 00. 

Let us begin with the band- width distribution Pb{w). 
We normalize the distribution by J PB(w)dw — 1 and 
/ °° wPb(w)cIw = 1. Each band width approaches zero 
at limit £ — > 00 [T5| . At a finite Fibonacci index 
£, we normalize the band widths by the mean of them 
and consider distribution Pjp (w) of the normalized band 
widths. Then we extrapolate Pb(w) at £ — * 00. 

Fig. 3 shows \n w-\n Pb {w) plots. As w — > 0, the 
distribution Pb{w) converges to zero. The results for 
different £ seems to be on the same line, which suggests 
the existence of band- width distribution at I — > 00. The 
distribution Pb(w) takes a form 

P B (w)~ w « (tu-»0), (13) 



with a positive real number a. This form remains to 
hold when / changes. The values of a for several cases 
are given in Table I. 

We next consider a region where w is larger than 0.5. 
In Fig. 4, we present distribution in w-hi Pb(w) plots. 
The linear behavior suggests that the distribution Pb{w) 
at large w is given by a form 

P B {w) - e~ f3w (w->oo), (14) 

with a positive real number (3. Estimated values of [3 for 
several cases are presented in Table I. We have changed 
the hopping coefficients t a and 4 for fixed /'s for several 
cases. The results indicate that a and (3 depend on t a 
and tb only through the value of /. 

Next we examine the ansatz that the distribution 
Pb(w) is described by a combined distribution P(w) — 
N a pw a e~P w with a normalization constant N a p for every 
value of w. This is the semi-Poisson form as was assumed 
for the Harper model The normalization conditions 
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2044!) 
144 
1225 
64 
22') 
_9_ 
16 
14400 



14161 
16 
9 

28224 



1.20(±0.03) 
1.28(±0.05) 
1.42(±0.07) 
1.71(±0.07) 
1.87(±0.13) 
1.91(±0.09) 
2.23(±0.17) 





0.38(±0.08) 
1.05(±0.08) 
1.07(±0.07) 
1.08(±0.06) 
1.18(±0.07) 
1.29(±0.08) 
1.50(±0.12) 



0.935(±0.004) 
0.824(±0.004) 
0.747(±0.004) 
0.689(±0.004) 
0.612(±0.002) 
0.555(±0.003) 
0.490(±0.002) 



P(w) 



TABLE I: Estimated values for critical indices a, (3 and S. 
Errors arise in fitting numerical data to the proposed distri- 
butions. 
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FIG. 5: Gap distribution lns-lnP G (s) for (t a ,t b ) - (1,2), 
/ = A and Fai = 17711. P G {s) ~ s" 689 . 



for the band-width distribution impose the relation be- 
tween the indices f3 = a + 1. Then the distribution re- 
duces to generalized form of the semi-Poisson distribution 
P(w) = { rt 1) ^ 1 w a e-( a+ V w . However this is different 
from the obtained behavior at w — > 0,oo. Thus we ex- 
clude the ansatz of the semi-Poisson distribution. 

We found that the gap distribution Pg(s) always di- 
verge for s — > 0. An example using the lns-lnPG(s) 
plot is shown in Fig. 5. We normalize the distribution 
as J Q Pg(s)cIs = 1. From Fig. 5, the gap distribution 
obeys an inverse power law 



Pg(s) 



(15) 



which diverges at the origin. Values of 6 for several cases 
are given in Table I. We have changed the hopping coeffi- 
cients t a and tb for fixed /'s for several cases. The results 
indicate that S depends on t a and tb only through I. 
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FIG. 6: Logarithmic normal distribution P(w) 

(log m -M) 2 
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with y, = 0, a = 0.6 



by splitting an interval into a definite number of subinter- 
vals with definite scales (see Appendix for details). We 
call the intervals of the multi-scale Cantor sets "bands" , 
while we call the removed intervals as "gaps" . 

Let us consider the band- width distribution of a multi- 



3 (Cantor) 



B 



w). As far as we know, 



Pb"™ (u) has not been appeared in the literature 



scale Cantor set, 

(Cantor) , 
B 

However, Kolmogorov's theory of turbulence |17j in- 
troduces the energy distribution measured at a defi- 
nite spatial scale, which is conceptually equivalent to 
P B CantOT \w) of our case. We derive P^ antOT \w) in the 
Appendix, and it turns out to follow a logarithmic normal 
distribution as a distribution in Kolmogorov's theory, 



^(Cantor) 
B 



(w) 



1 



N(o,na 2 )(\nw), (n-s-oo), (16) 



where N, 



(,:) 



~? r ~ and <t 2 is the variance 



of the logarithm of the scales to generate a multi-scale 
Cantor set and n is the number of steps to produce the 
set. Thus for a multi-scale Cantor set, P^ antm \w) be- 
comes broader for large n and does not converge to a 
normalizable distribution. This contrasts with the be- 
havior of Pb{w) for the Fibonacci model. 

As shown in the Appendix, we can modify the con- 
struction of a multi-scale Cantor set so that p^ Cantor ) ^) 
has a constant variance by including a finite size correc- 
tion to the scales at each splitting step. When we write 
the width of a band produced at step k as Wk = e~ hek , 
ek would receive a correction proportional to ^d*-. Then 



^(Cantor) 
B 



(w) = —N^^Qxiw). 



(17) 



IV. COMPARISON WITH THE MULTI-SCALE 
CANTOR SETS 

We compare our results in the previous section with 
properties of multi-scale Cantor sets which are produced 



A graph for a logarithmic normal distribution is shown 
in Fig. 6. At the origin, the logarithmic normal distri- 
bution approaches zero faster than power law. Also at 
the tail, it does not decay at an exponential rate. Thus, 



the behavior of P 



(Cantor) 
B 1 



w) of this multi-scale Cantor 
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set with the finite size corrections shows differences from 
those of the Fibonacci model. 

These differences may be rooted in the degree of self- 
similarity. The multi-scale Cantor sets are genuinely self- 
similar : the same splitting procedure is applied to every 
interval at each step. This leads to a logarithm nor- 
mal contribution. On the other hand, the splitting of 
the energy levels of the Fibonacci model is different for 
each band and it has different scales for different bands 
[Til ITU . In , it has been observed that numerical val- 
ues of scales of the band's center coincide with eigenval- 
ues of the linearized KKT map at the hyperbolic points of 
the corresponding periodic orbit and the splitting of the 
bands is explained by the Smale's horseshoe structure. 
This observation can be applied to other periodic orbits 
corresponding to other bands. Then, at least, some sub- 
sets of the bands are self-similar and should satisfy the 
band-width distribution like the multi-scale Cantor sets. 
Still a finite number of periodic orbits only lead to the 
sum of logarithmic normal distributions and cannot ex- 
plain fully the behavior of the Fibonacci model. 

Then our comparison implies that, to understand the 
behavior of Pb(w) of the Fibonacci model by the lan- 
guage of dynamical system, it is necessary to include fi- 
nite size corrections as well as an infinite number of pe- 
riodic orbits. Such a fine structure of fractal is not taken 
into account in multifractal analysis of 0, which char- 
acterizes fractal objects by the statistical distribution of 
scales. As proved in 0, the energy spectra of the Fi- 
bonacci model belong to a class of a general Cantor set. 
Thus they belong to a different class of Cantor sets to 
the multi-scale Cantor sets. 

We also compare the results for the gap distribution. 
The gap distribution of the multi-scale Cantor set follows 
an inverse power law near the origin s — > 



P 



(Cantor) 



G 



00 



-6 



S = D 



H 



1. 



(18) 



where Dh is the Hausdorff dimension. This is because 
the length of gaps of a multi-scale Cantor set is directly 
related to the length of the removed elements since the 
bands and gaps are complementary in a fixed interval. 
Although the Fibonacci model have an inverse power 
law like Ea. l|18|) for the gap distribution, the relation 
8 = Dh + 1 does not hold for the Fibonacci model be- 
cause of the value of the Hausdorff dimension ~ 0.685 
[l4| . This discrepancy is due to the splitting of the en- 
ergy levels as observed in Fig. 1: the energy levels of 
the Fibonacci model have fluctuating endpoints, and the 
bands and the gaps are not complementary of a fixed in- 
terval, and hence have different scaling properties. On 
the other hand, for the Harper model, the numerical re- 
sults of 0, 0] seem to sustain the relation 6 = Dh + 1 • 
From the same reasoning, there seems to exist no expla- 
nation for it actually. 



V. CONCLUSION 

Spectral properties of the Fibonacci model are ex- 
amined numerically by using the Kohmoto-Kadanoff- 
Tang map. We investigated the band-width distribution 
Pb(w) and the gap distribution Pg(s). We find a power 
law Pb{w) ~ w a near the origin and an exponential de- 
cay Pb(w) ~ e~P w toward infinity (a,f3 > 0). Also P G (s) 
diverges with an inverse power law s~ s near the origin 
(5 > 0). The indices turn out to depend on the cou- 
pling constants through the constant of motion in the 
Kohmoto-Kadanoff-Tang map. 

We compare the band- width distribution Pb(w) and 
the gap distribution Pg(s) of the Fibonacci model with 
those of multi-scale Cantor sets. As the distribution 
of the "band" widths P^ Cantor - ) (w) for multi-scale Can- 
tor sets always obeys a logarithmic normal distribution, 
the behavior of Pb(w) is rather different from that of 



^(Cantor) 
B 1 

rj(Cantor) | 
G 



w). Also, while the distribution of the "gaps" 
s) of the multi-scale Cantor sets is an inverse 
power law s~ s near the origin (5 > 0) and there is the re- 
lation between the index S and its Hausdorff dimension, 
the relation does not hold for the Fibonacci model. Thus 
the critical level statistics is useful to quantify differences 
between the energy spectra of the Fibonacci models and 
multi-scale Cantor sets which are hidden in the language 
of multifractal analysis. 
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APPENDIX: DISTRIBUTIONS FOR THE 
MULTI-SCALE CANTOR SETS 

The general definition of a Cantor set is a closed set 
with no isolated points and whose complement is dense. 
Here we restrict ourselves to the so-called multiple-scale 
Cantor sets. A multi-scale Cantor set is a set which is 
obtained by the following procedure : we transform the 
interval / = [0,1] into a finite number of subintervals 
with definite scales and apply the same transformation 
to subintervals, and the limit set is called a multi-scale 
Cantor set. 

A standard example is the classical Cantor set : We 
split an interval into three equal pieces and remove the 
middle one. Two remaining intervals are [0, ^] and [|, 1]. 
One can continue the same transformation for these in- 
tervals. The number of intervals which we call "band" 
at step n is given by N — 2™ with length W = 
The removed intervals, which we call "gap" , have lengths 
S = §, 52, • • • , .^r- As n goes to oo, we end up with an 
infinite number of infinitesimal intervals. The Hausdorff 
dimension of the set is given by Dh = = jf§ • 
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For a multi-scale Cantor set, let n be the number of 
steps to generate the set and N be the number of bands : 
N ~ a n with a number of split intervals a. Let Wi, Si (i = 
1, • • • , N) be the length of bands and gaps respectively. 
We normalize Wi so that their mean is equal to 1. The 
mean of Wi at step n is given by W n = jfY^iLi^i- The 
normalized band-widths are Wi 



w„ 



For the Cantor set with one scale described above, N 



2™,Wi = 3" 
1, we have 



N W = Eili Wi='g,W = -±; and w t = 



2H 

3" 



P B Caator \w)=S(w-l). 



(19) 

This means that there is no distribution of scaling vari- 
ables and is readily generalized to a scale other than 
■|. On the other hand, the gap lengths are Si = 
Ii^)'*')^- The number of the gaps of ^ (k = 
1, • • • , n) is 2 k . Thus the distribution of the gaps is given 
by 



p (,Cantor) (s) 



(20) 

The power is related to the Hausdorff dimension because 
the lengths of the gaps consist of those of the band- 
widths. 

Let us next consider a two-scale Cantor set which has 
two scales r\ , r% > with r± + ri = po < 1 . Namely 
we split an interval into three intervals with length r\ , r 2 
and t = 1 — po, and remove the interval with the length t. 
Then we continue the same procedure to the remaining 
intervals. We put N — 2", p = p a /2 and define x = 

x/p,y = y/p- 

For the bands, we have 



w. — r n r *»-l_ r n-2 2 r ,r n ^ 1 r n 

yy% — I lil l 1 2,' i '2> >'1'2 >'2: 



(21) 

where each of them has ( ra ) terms. Then, the mean of 
the band- widths is given by 



TV A = 



(ri+r a )» = pj, 



A = p™. 
The normalized band-widths are 



(22) 



Since the binomial distribution (™) converges to a normal 
distribution as n — > oo, when n is taken as a discrete 
time, a; can be regarded as a discrete Brownian motion 
(a binomial model of it). To get P B Cantol \w) for n — > oo, 
let us recall the central limit theorem: 
[Theorem] (The Central Limit Theorem) 
Let Xk be independent stochastic variables following 
an identical probabilistic distribution. Let E(X\) = 



p, VariX^ 



lim Prob 



X, 



Then 
+ X 2 - 



x„ 



np 



< x 



awn 



N 0< i(y)dy, N^ a 2 } (y) 
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In our case, we define Xk 
Xk take values 

"±d,Var{X k ) = 
applies : 



In 



(Wo 



1) 



u = lnri,d = lnr 2 , and E(Xf.) = 
, Then the central limit theorem 



lim Prob 



In 2^ 

111 W 



n{u+d) 
2 



j(u — d)^/n 



< X 



(24) 



This 

iV, 



means 



that In -pj^ obeys the distribution 
^ . Returning to the normalized variables 

-. The band-width 



(f(u+d) 

tDj, we similarly define x k = In 
distribution is obtained as 



^(Cantor) 
B 



W (0: 



j(lnw), (n — > oo). 

This is a logarithmic normal distribution whose variance 
increases linearly with n. 

The distribution (|25|l can be generalized to Cantor sets 
with arbitrary number of scales as long as the mean and 
the variance of the scales are finite. In this case, the cen- 
tral limit theorem applies. P B Caatm \w) always obeys a 
logarithmic normal distribution with a variance propor- 
tional to n : 

1 



P, 



(Cantor) 



B 



(w) N(p tTUT z)Qjiw), (n->-oo), (26) 



where a 2 is the variance of the logarithm of the scales. 

(w) with a definite variance, one 



(Cantor) 



To obtain P B 

can modify the transformation of splitting a interval. For 
simplicity, we describe it for the two-scale Cantor set. We 
divide the interval [0, 1] into three intervals with lengths 
ri, r 2 and 1 — p with n, r% > 0, ri + r<z = po < 1, and 
remove the interval with the length 1 — po- We continue 
the procedure for the remaining intervals but, at step k, 



we always replace r% and r 2 with r%(k) = r± 



and 



•/k 



When we write the width of a band 



r 2 (k) = r_ 

produced at step k as Wk — e~ fcCfc by an exponent e^, 
the factor is amount to change et by a finite size correc- 
tion Then the distribution becomes a logarithmic 
normal distribution with a constant variance: 



P 



(Cantor) 



B 



w C ' — — 



(law). 



(27) 



This can be readily generalized to the cases with more 
than two scales, leading to the distribution 



P. 



(Cantor) 



l 



B 



(to) = -% ff 2)(lnM)), 



(28) 



where a 2 is the variance of the logarithm of the scales. 

For the gap distribution of the multi-scale Cantor sets, 
the formula p^ Cantc ") (s')ds' ~ s~ Dh with their Haus- 
dorff dimension Dh holds generally because the lengths 
of the gaps are directly connected with those of the band- 
widths. Thus the gap distribution for the multi-scale 
Cantor sets is given by 

p (Cantor) (s)rvjs _D„-l (2Q) 
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